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The emergence of an effective field theory out of equilibrium is studied in the case in which a 
light field -the system- interacts with very heavy fields in a finite temperature bath. We obtain 
the reduced density matrix for the light field, its time evolution is determined by an effective action 
that includes the influence action from correlations of the heavy degrees of freedom. The non¬ 
equilibrium effective field theory yields a Langevin equation of motion for the light field in terms 
of dissipative and noise kernels that obey a generalized fluctuation dissipation relation. These are 
completely determined by the spectral density of the bath which is analyzed in detail for several 
cases. At T = 0 we elucidate the effect of thresholds in the renormalization aspects and the 
asymptotic emergence of a local effective field theory with unitary time evolution. At T ^ 0 new 
“anomalous” thresholds arise, in particular the decay of the environmental heavy fields into the 
light field leads to dissipative dynamics of the light field. Even when the heavy bath particles are 
thermally suppressed this dissipative contribution leads to the thermalization of the light field which 
is confirmed by a quantum kinetics analysis. We obtain the quantum master equation and show 
explicitly that its solution in the field basis is precisely the influence action that determines the 
effective non-equilibrium field theory. The Lindblad form of the quantum master equation features 
time dependent dissipative coefficients. Their time dependence is crucial to extract renormalization 
effects at asymptotically long time. The dynamics from the quantum master equation is in complete 
agreement with that of the effective action, Langevin dynamics and quantum kinetics, thus providing 
a unified framework to effective field theory out of equilibrium. 


I. INTRODUCTION 

Effective field theory is a powerful organizational principle to describe phenomena below some energy scale, or 
alternatively on large spatio-temporal scales, and is ubiquitous across fields. Several applications of effective field 
theory have become the pillars of fundamentally important paradigms, for example: universality in critical phenomena 
emerges at long wavelengths after coarse graining over short wavelength degrees of freedom a la Wilson: the Landau- 
Ginsburg theory of phase transitions classifies universality classes in terms of few coarse grained order parameters and 
their symmetries [lj, the Landau-Ginsburg theory of superconductivity emerges after integrating out the fermionic 
quasiparticles leading to an effective theory of the superconducting order parameter [2], and hydrodynamics is a 
description of long wavelength collective flow that emerges after coarse graining over small scales and is valid on 
scales much larger than a microscopic mean free path. These are but a few historically and conceptually important 
examples of effective field theories. In particle physics[^-@| effective field theory provides a systematic characterization 
of phenomena below some energy scale and is an important tool to incorporate the physics associated with the degrees 
of freedom beyond such energy scale in a consistent and systematic manner. These concepts have been extended to the 
realm of early Universe cosmology to describe cosmological perturbations [2, H|, large scale structure formationjsj and 
could ultimately underlie the description of inflationary cosmology in terms of an inflaton scalar field as an effective 
dynamical degree of freedom well below the Planck scale [ic|. The usual approach to effective field theory begins by 
recognizing the operators that could enter in an effective Lagrangian based on internal and space-time symmetries 
in a suitable derivative expansion, where higher derivatives and higher dimensional operators are suppressed by 
inverse powers of the (high) energy scale[sj-Q, such a description is manifestly local and yields unitary time evolution 
of observables. However, it is not always the case that integrating over “hard” scales the resulting low energy 
effective field theory is local, an important counterexample is the “hard thermal loop” effective field theory at finite 
temperatureflU PHI. At finite temperature, Landau damping is a medium process that yields spectral densities in loop 
diagrams that feature support below the light cone and leads to non-local terms in the effective action[ll| and long 
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time tails in the dynamics [l^. Fundamentally, effective field theory must be understood as emerging from “tracing 
over”, “integrating out” or “coarse graining” high energy degrees of freedom with fast dynamics and short wavelength 
fluctuations, and it describes the influence of the high energy degrees of freedom over the low energy (slow, long 
wavelength) degrees of freedom. Cast in this manner, effective field theory is another manifestation of the quantum 
open system approach to studying the effect of an environment upon the dynamics of a system [l4l - [l7l ] pioneered with 
the study of quantum Brownian motion [l8l - {20| . This approach begins by considering the time evolution of the full 
density matrix of the system coupled to the environment and tracing over the environmental (or “bath”) degrees of 
freedom leading to a reduced density matrix for the system, whose time evolution includes the influence action from the 
environment |18( upon the system. Within this context the effective action arises as the sum of the system’s action and 
the influence action. The effective equations of motion for the system variables become a quantum Langevin equation 
with a dissipative term and stochastic force that obey a generalized fluctuation dissipation relation [2l| - l23j . Quantum 
Brownian motion of a particle coupled linearly and non-linearly to general environments has been studied thoroughly in 
refs.[24[ and a comprehensive and in-depth discussion of non-equilibrium phenomena and in particular generalizations 
of quantum Brownian motion is available in refs. [25l - (27i |. Decoherence and effective stochastic dynamics emerging 
from tracing over short wavelength degrees of freedom are also of fundamental importance in cosmology [2^-0]. 

The influence action approach has recently been argued to provide an effective description of the dynamics of 
long-wavelength fluctuations when combined with a Wilsonian approach to coarse graining the short wavelength 
components[35J which are taken as an environment or bath. An alternative approach to quantum open systems relies 
on the quantum master equation [3-[I3| for the reduced density matrix. This approach has recently been advocated 
in cosmologyj36] under the assumption that the environment only features short time (delta function) correlations. 

Goals : In this article we study the emergence of an effective field theory description out of equilibrium in the case 
in which a light field -the system- interacts with very heavy fields, taken to be the “environment” and integrated 
out or traced over. In particular we focus on analyzing in detail the influence of correlations of the heavy degrees 
of freedom, upon the non-equilibrium dynamics of the light degrees of freedom. We seek to elucidate in a direct 
manner the relationship between the influence action, stochastic and quantum master equation approaches without 
assumptions on environmental correlations, focusing precisely on how the spectral properties of the environment lead 
to the different dynamical time scales of the low energy degrees of freedom. We address the following questions: (I): 
In quantum field theory there are thresholds to excitation of heavy degrees of freedom, a local and unitary effective 
action emerges from integrating out heavy (or high energy) degrees of freedom when the energy and momentum 
associated with the dynamics of the light fields is well below this threshold. How are these threshold effects manifest 
in the influence action, stochastic and quantum master equation description?. (II): When the heavy degrees of 
freedom form a thermalized plasma through their mutual interaction, there are new “in medium” corrections to the 
spectral density of correlators of the heavy degrees of freedom that may lead to dissipative processes in the influence 
action of the light degrees of freedom even when the T = 0 thresholds correspond to high energy. This is the case 
in “hard thermal loops” HE HI as a consequence of Landau damping. Do these new, “in medium” contributions 
to the spectral density lead to dissipative dynamics of the light field?. (Ill): does the light field thermalize with 
the heavy coarse grained degrees of freedom?. If the light field tliermalizes with the “bath” of heavy particles, do 
the influence functional and quantum master equation approaches agree with the dynamics of thermalization from 
quantum kinetics?. (IV): Heavy environmental fields may actually decay into the lighter species considered as the 
“system” when they are coupled, how does the decay of the environmental degrees of freedom affect the dynamics 
of the light fields?. (V): How do renormalization aspects emerge from the non-equilibrium descriptions: influence 
action, stochastic and quantum master equation approaches?. 

Summary of results: To answer these questions we focus on the case of a light scalar field <fi of mass coupled 
to either just one scalar heavy field \ with M or two heavy scalar fields XuX 2 with a hierarchy of masses 

M i > M 2 77i0. The heavy fields generically denoted as x are the environmental degrees of freedom, which are 
integrated out or “traced over”, the light field <j> constitutes the “system”. We study couplings of the form 0[\j] 
where J[(j)\ and O[xj ] are generic polynomials in the respective fields. We focus in detail on two types of couplings 
from which we draw more general conclusions: (i): g<f> 2 Xi (ii) : gf>X 1X2 both at T = 0 and T / 0. The first type 
yields an effective “current-current” interaction in the low energy limit and illustrates how this familiar local effective 
field theory emerges in the influence action and quantum master equation approaches in the long time limit. The 
second type introduces a rich spectral density with “anomalous” thresholds to excitation of the degrees of freedom of 
the environment as a consequence of in medium effects. As a consequence of the heavy field with M\ decaying into 
the light field <j> this spectral density features support on the mass shell of the light particles leading to dissipative 
phenomena with a wealth of dynamical scales. The main results are: 

• We obtain the time evolution of the reduced density matrix and the non-equilibrium effective action by tracing 
over the “bath” degrees of freedom to order g 2 in terms of the various correlation functions of the “bath”, 
which are determined by the spectral density in all cases. For the case g(f 2 x (single heavy field), we analyze 
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the transient dynamics and obtain the “current-current” effective theory in the long time limit. In the case of a 
thermal bath of two heavy fields XuX 2 with Mi > M 2 we find the (one loop) spectral density for T / 0, 

it features several thresholds below the two particle threshold at T = 0, one describes Landau damping and the 
other describes the decay of the heavier bath field into the light fields. These two “anomalous” thresholds lead 
to a dissipative contribution to the <f> effective action, in particular we find that the decay Xi —> X 2 4> leads to a 
purely dissipative contribution to the effective action which describes the thermalization of the light field with 
the bath. 

• From the non-equilibrium influence action we obtain a semiclassical stochastic description for the light field, 
the equation of motion for the non-equilibrium average is a generalized Langevin equation with a stochastic 
Gaussian noise term and a non-local dissipative kernel- hence <f> is a light Brownian quantum field-, the noise 
and dissipative kernels obey a generalized fluctuation-dissipation relation. For the coupling g<j>x 1 X 2 the noise is 
additive but with a colored spectrum that reflects the underlying temporal correlations of the bath, for more 
general couplings g J\<j\x 1 X 2 with J[<fi\ ^ <f> we find multiplicative noise. For the case g<f>x 1X2 we obtain the 
solution of the Langevin equation and obtain the correlation functions of the light field, for T/Owe find that 
the decay xi X 2 </ ) leads to the thermalization of <j) and obtain the time scale for thermalization. We show 
that the dynamics of thermalization obtained from the effective action is in complete agreement with a quantum 
kinetic description of <f> thermalization. As T —> 0 only the high energy two particle threshold remains in the 
spectral density of the bath, which is responsible for renormalization effects such as wave function and mass 
renormalization. In this case a local and unitary effective field theory emerges in the asymptotic long time 
limit. 

• We obtain the quantum master equation up to 0(g 2 ) and show that its solution in the field basis is precisely 
the influence function obtained by tracing over the bath. In the case of linear coupling to the bath we obtain 
the Lindblad form of the master equation under clearly specified approximations. The Lindblad form contains 
a Hamiltonian and a dissipative term, both feature time dependent coefficient functions. The time dependence 
of the dissipative functions is crucial to understand renormalization effects. The quantum master equation 
describes thermalization in complete agreement with the influence function, the Langevin stochastic description 
and the quantum kinetic equation. Keeping the time evolution of the dissipative coefficients in the Lindblad 
form allows to extract the dynamics of the “dressing” of the bare states and wave function renormalization in the 
asymptotic long time limit, again in complete agreement with the influence function and Langevin descriptions 
even in the T = 0 case when only the high energy thresholds are present in the spectral density of the bath. 


II. THE NON-EQUILIBRIUM EFFECTIVE ACTION 

We study the non equilibrium effective action for the case of interacting scalar fields, spinor and vector fields may be 
included, while technically more involved their treatment follows without major conceptual difficulties. We consider 
a bosonic field <t> of mass referred to as the “system” and either a single massive scalar field x of mass M x 
or two scalar fields collectively denoted by Xa ; * = 1,2 with the hierarchy of masses Mi > M 2 and a generic 

interaction between <fi and Xi- The Lagrangian density is given by 

£[<!>{x),Xa{x)] = c 0 ^[<l>(x)\ + C 0 , x \xa{x)\ -g jy>(x)\ 0[xa(x)] (2.1) 

where the “current” J[4>\ and O[xo] are in general non-linear function(als) of the respective fields, and Co,fax are 
the free field Lagrangian densities for the fields <j> and Xa , again self-interactions may be included at the expense of 
technical complications but with no conceptual difficulties. 

Specifically we consider the following cases: 

• Only one field x and J[<f>(:x)\ = (fr(x), namely with the interaction g (jr(x) xi x )- This interaction describes the 
exchange of a massive “vector boson”, x whereas «/[(/>] = <f> 2 describes a bilinear “current”, as in the coupling 
between gauge bosons and fermionic degrees of freedom. This model allows us to understand the emergence 
of a “local Fermi” theory in the limit where the frequency and momentum transferred by the current is much 
smaller than the mass of the (vector) boson. 

• Two different fields Xi ,2 with Mi > M 2 and J\<p(x)\ = 4>{x), namely with interaction g Xi(x) X 2 (x). 

This case allows us to obtain a Langevin equation of motion for the system field (j> describing its non-equilibrium 
dynamics as a “Brownian” field. This case will also lead to a detailed understanding of thermalization and 
dissipative processes by interactions with heavy fields. 
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• Two different fields with the same hierarchy as the previous case but now with J\<j){x)\ = (jr{x) and interaction 
g cf> 2 (x) Xi(x) X 2 {x). This case will showcase important renormalization aspects and highlights the limitations 
of a local description when the “bath” or environmental fields xi ,2 form a plasma. 


While we are ultimately interested in obtaining an effective quantum field theory by tracing out “heavy degrees of 
freedom” in cosmology, in this study we focus on Minkowski space time and consider that the fields Xi are treated as 
a bath in equilibrium assuming that the bath fields are sufficiently strongly coupled so as to guarantee their thermal 
equilibration. These fields will be “integrated out” yielding a reduced density matrix for the field <j> in terms of an 
effective real-time functional, known as the influence functional[18] in the theory of quantum brownian motion. The 
reduced density matrix can be represented by a path integral in terms of the non-equilibrium effective action that 
includes the influence functional. This method has been used previously to study quantum brownian motion 
and for studies of quantum kinetics beyond the Boltzmann equation }37 h39| . 

Consider the initial density matrix at a time t = 0 to be of the form 


P{ 0) = /5fl>(0)<g>p x (0). 


( 2 . 2 ) 


The initial density matrix of the Xo fields will be taken to describe a statistical ensemble in thermal equilibrium at 
a temperature T = 1//3, namely 


Px( 0) = e~ f>H °*, (2.3) 

where Hq x {xo) is the free held Hamiltonian for the fields Xa- We will now refer collectively to the set of fields Xa 
simply as x to avoid cluttering of indices. 

The factorization of the initial density matrix is an assumption often explicitly or implicitly made in the literature, 
it can be relaxed by including initial correlations, we will not consider here this important case, relegating it to future 
study. The initial factorization entails that the time evolution of the system can be described as “switching-on” the 
coupling between fields at the initial time. This will result in transient dynamics, however, we will focus on the long 
time evolution. 

In the held basis the matrix elements of p c / > ( 0) are given by 

= PtAfatf) ; (x|Px(°)lx') = Px,o(x;x'), ( 2 - 4 ) 

we emphasize that this is a functional density matrix as the held has spatial arguments. The density matrix for <f> 
represents an initial out of equilibrium state or ensemble. 

The physical situation described by this initial initial density matrix is that of a held (or helds) in thermal equilibrium 
at a temperature T = 1//3, namely a heat bath, which is put in contact with another system, here represented by 
the held <j>. Once the system and bath are put in contact their mutual interaction will evolve the initial state out of 
equilibrium because the initial density matrix does not commute with the total Hamiltonian. 

To obtain the effective quantum held theory out of equilibrium for the light held <f> we will evolve the initial density 
matrix in time and trace over the “bath” degrees of freedom, leading to a reduced density matrix for <f >. Once we 
obtain the reduced density matrix for the held </> we can compute expectation values or correlation functions of this 
held. 

The time evolution of the initial density matrix is given by 


p(f) = U{t)mU~ l (t) 1 (2.5) 

where 

U{t) = e ~ lHt . 

The total Hamiltonian H is given by 

H = Hoc/, + H 0x + i?/(0, x) 5 #/(<£, X) = 9 j d 3 x J[(j>{x)} 0\x{x )\, 
and Hod,, Hq x are the free held Hamiltonians for the respective helds. 

= <^>/; X/|C/(t)p(O)J7 _1 (t)|0^; X/) 

= J DfaDxiD^Dx'i (<ff;xf\u(t)\(j>i]Xi)PtAfa&i)Px,o(xi;x'i) <^;Xil^ -1 (*)I‘A/;X/) ( 2 -8) 


( 2 . 6 ) 

(2.7) 
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The f D(f> etc, are functional integrals where the spatial argument has been suppressed. The matrix elements of the 
time evolution forward and backward can be written as path integrals, namely 


(fa',Xf\U(t)\fa,Xi) = J 

vfav x + e iidixC[<t,+ ' x+] 

(2.9) 

(fafaifa^fafaflx'f) = J 

VfaVx~e~ i fa 3xC ^~' x ~ ] 

(2.10) 

where we use the shorthand notation 



S* zs L 

dt J d 3 x , 

(2.11) 

C[4>, x] is given by (12.11) and the boundary conditions on the path integrals are 


fa(x,t = 0) = fax) ; 

fa(x,t) = fif(x), 


X + (x,t = 0) = Xi(x) ; 

X + (x,t) = Xf(x), 

(2.12) 

<t> (x,t = 0) = fax) ; <j> (x,t) = <t>' f (x), 


x (x,t = 0) = fax) ; 

X (x,t) = x'f(x)- 

(2.13) 


(2.14) 

The field variables fa,"fa along the forward (+) and backward (—) evolution branches are recognized as those 
necessary for the Schwinger-Keldyshpjl [25], [|o|, [IjJ closed time path approach to the time evolution of a density 
matrix. 


A. Tracing over the “bath” degrees of freedom: reduced density matrix 

The reduced density matrix for the light field <j> is obtained by tracing over the bath (xi) variables, namely 

P r (<t>f,<t>'f\t) = J Dxf , (2.15) 

we find 

P r ( ( t ) ffa'f'fl) = J D 4>i D 4>'ifa^ffa'f\fafaAp4>{fafa 0), (2.16) 

where the time evolution kernel is given by 

T[4> f , fa fa, fat] = Jvfa JVfa e i:F[J+ ’ J ~ ] (2.17) 

with the following boundary conditions on the forward (fa) and backward (fa) path integrals 

fa (x , t = 0) = <j>i (x) ; fa (. x , t) = </>/ (x) 

fa(x,t = 0) = fax) ; fa(x,t) = fa(x). (2.18) 

J-[J + \ J - ] is the influence action where fa = J[fa] given by 

fal J+ ; J ~] = j Dx . J Dx'faxf J T> X + J V X ~ e ifdix { Colx+] - 9Ji * +] ° ix+] } e -i/^*^o[x-]- 9 J[0-]O[ x -]] Px ( XijX >, 0 ) 

(2- 19 ) 

the boundary conditions on the path integrals are 

X + (x,t = 0) = Xi(x) ; X + (x,t) = Xf(x) ; X~(x,t = 0) = x'(f) ; X~(x,t) = x!f(x) ■ (2.20) 

In the above path integral defining the influence action, ./[(Z^] acts as an external source coupled to the composite 
operator 0(x), therefore, it is straightforward to conclude that 


fafa’J 1 =Tr \u(t;fa)p x (0)U- 1 (t;J-) 


( 2 . 21 ) 
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where J ± = and U(t\ J±) is the time evolution operator in the x sector in presence of external sources J ± 

namely 


U(t-J + ) = T(e~ i lo H x[ J+ (.t')W^ . u~ 1 (t;J-) = T(e i ^ H xlJ = H 0x + J d 3 *J[^ ± (t)] 0 (x) 

( 2 . 22 ) 

and T is the anti-time evolution operator as befits evolution backward in time. The calculation of the influence action 
is facilitated by passing to the interaction picture for the Hamiltonian H x [J(t)], defining 

U{t;J ± ) = e- iH °x t U ip (t;J ± ) (2.23) 

and the e ±lH ° x 1 cancel out in the trace in (12.211) . Now the trace can be obtained systematically in perturbation theory 
in g. Up to 0(g 2 ) we find 


T[J + ,J ] = — g J d 4 x(^J + [x\ — J [x]^ (0(x)} + 

J d 4 x i J d 4 x 2 < J + [x\\ J + [x 2 ] G+ + (xi - x 2 ) + J _ [xi] J~[x 2 \ G~~(xi - x 2 ) 


- J + [xi]J [x 2 ]G+ (xi — x 2 ) — J [xi] J + [x 2 ] G c + (xi - x 2 )|. (2.24) 

In this expression J^x] = J[<^ ± (x)], and the connected correlation functions are given by 

G-+(X! - x 2 ) - {0(x 1 )0(x 2 )) - <0(x 1 ))(0(x 2 )> = G>(X! - x 2 ), (2.25) 

G+-(xr - x 2 ) = (0(x 2 )0(xi)) - <0(x 2 ))<0(x 1 )) = G<(x! - x 2 ), (2.26) 

G+ + (xi - x 2 ) = G>(xi - x 2 )0(ti - t 2 ) + G<(xi - x 2 )0(t 2 - ti), (2.27) 

G c (xi — x 2 ) = G^(xi — x 2 )0(t 2 — ti) + Gf(x\ — x 2 )0(ti — t 2 ), (2.28) 

in terms of interaction picture fields, where 

((''')) = Tr(- • • )Px(0) • (2-29) 

Furthermore, for the case of hermitian operators O as considered here it follows that 

G^ (xi — x 2 ) = Gf(x 2 — xi). (2.30) 

The effective action out of equilibrium is given by 

Seff[4> + ,<l>-] =j*dt J d 3 x\cv[(t> + }-Cv[r}\+nJ[4> + },J[r]]- (2.31) 


Since p x (0) = e _/3i4 ° x the expectation value (0(x i)) is independent of space-time coordinates. Therefore the first 
term on the right hand side of (12.241) can be cancelled by normal ordering the composite operator 0(x) in the initial 
density matrix of the x field, namely 


0(x) -> 0(x) - (0(x)) • (2.32) 

This normal ordered operator features vanishing expectation value in the initial density matrix. This is tantamount 
to introducing a counterterm in the Lagrangian density that cancels the first term on the right hand side of (12.241) 
and the terms (G(xi j2 )) in the connected correlation functions (12.25112. 251) . hence in what follows we suppress the 
subscript “c” in the correlation functions. 

With the purpose of comparing the influence functional approach to the quantum master equation developed in 
section m and to exploit the spectral representation of the correlation functions, we re-write the influence action 
(12.241) solely in terms of the two correlation functions G^. This is achieved by implementing the following steps: 

• In the term with J + (xi) J + {x 2 ): in the contribution G < (xi — x 2 )0(t 2 — t\) (see eqn. (12.271) 1 relabel t\ -H- t 2 
and use the property (12.301) . 
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• In the term with J {x\)J (x 2 ): in the contribution G > (x 1 — x 2 )0(f 2 — G) (see eqn. (12.2811 ) relabel G -H- t 2 
and use the property (12.301) . 

• In the term with J + {x\) J~ (x 2 ): multiply G < (x 1 — X 2 ) by 0(G — £ 2 ) + 0(i 2 — G) = 1 and in the term with 
0(i 2 — G) relabel G t 2 and use the property (12.301) . 

• In the term with J~(x\) J + (x 2 ): multiply G > (xi — x 2 ) by 0(G — t 2 ) + 0(f 2 — G) = 1 and in the term with 
0(t 2 — G) relabel G «-»■ t 2 and use the property (12.301) . 


We find 


T[J + ,J ] = Wj d 3 x 1 d 3 x 2 J dti J dt 2 | J + (x 1 ,ti)J + (x 2 ,t 2 ) G > (xi - x 2 ) + J (xi,ti)J (x 2 ,t 2 ) G < (x 1 - x 2 ) 


- J+(xi,ti)J (x 2 ,f 2 ) G < (x\ — x 2 ) — J (x 1 ,ti)J + (x 2 ,t 2 )G > (x 1 - x 2 ) 


(2.33) 


where G§ are given by eqns. (12.2512.261) . This is the general form of the influence function up to second order in the 
system-environment coupling. 

We note that 


d_ 

dt 


p iF[J + ,J-} l _ 


= ~9' 


J d 3 x\d 3 x 2 J dt 2 | J + (xi,t)J + (x 2 ,t 2 ) G > (fi - x 2 ,t - t 2 ) + 


J (xi,t)J (x 2 ,i 2 ) G < (x 1 - x 2 ,t - t 2 ) - J + (fi,t) J (x 2 ,t 2 ) G<(xi - x 2 ,t - t 2 ) 
J~(xi,t)J + (x 2 ,t 2 ) G > (x l - x 2 ,t - t 2 ) 1 J e ^[J + ,J ll (2.34) 


this result will be of paramount importance when comparing to the quantum master equation in section (|Vl) . 


B. Spectral densities of the environment: 

The dynamics and dissipative processes depend on the correlation functions of the environment and crucially on their 
spectral density. In appendix (A) we show that the correlation functions (12.25112.26]) (and consequently (12.2712.281) 1 
can be written in terms of a spectral representation, namely 

G% -x')= [ p^(k 0 ,k) e - ifc o(t-t') e iHx-S') (2.35) 

J (2tt) 4 


where 


p > (k 0 ,k) = p(k 0 ,k)[l + n(k 0 )] ; p < (k 0 , k) = p(k 0 , k)n(k 0 ) ; n(k 0 ) = ^_ 1 ; P = 1 /T (2.36) 

where p(k 0 ,k) is the spectral density and T is the temperature of the bath described by the y field(s). We analyze 
the different cases: (a): one single y field of mass M x m^, and 0(x) = X, and (b): two different fields yi,y 2 with 
0(x) = Xi X 2 and masses M 1 > M 2 m^. 

The spectral densities for the different cases are obtained in detail in appendix m, here we summarize their main 
features for the different cases. 

Case a): One y field, 0(x) = X 

In this case it is straightforward to find (see Appendix B) 


p{ko,k ) 


7r 

w k 


S(k 0 - wk) - S(ko + Wk) 


where M is the mass of the x held. 

Case b): Two y fields, G(y) = X1X2 


Wk = \Jk 2 + M 2 , 


(2.37) 





















In Appendix B we provide the derivation of the spectral density at finite temperature T for this case and show that 
it is of the form 


p(ko, k; T) = p LD {k 0 , k; T) LD 0(-Q 2 ) + p D (k 0 , fc; T) ©((Mi - M 2 ) 2 - Q 2 ) 0(Q 2 ) 

+ p 2x (fc o ,fc;T)0(Q 2 -(M 1 + M 2 ) 2 ) ; Q 2 = k 2 0 - k 2 (2.38) 

The expressions for the different p's are given explicitly in appendix B and merit discussion. 

• The contribution pLD(ko, k; T ) with support below the light cone ( Q 2 < 0) corresponds to the process of Landau 
damping, this is a medium dependent contribution that vanishes in the T —> 0 limit. It describes collisionless 
damping in a medium as a consequence of dephasing. While this contribution vanishes on the mass shell of the 
light field </> at Q 2 = ro|, it does contribute in the long time, long wavelength limit, as discussed in referencefH] 
it is responsible for power law long time tails. 


• The contribution /9_o(fco, k; T ) also describes a process solely available in the medium corresponding to the decay 
Xi —•► X'20 (since Mi > M 2 m^), this interpretation will be discussed in detail below. This part of the spectral 
density also vanishes for T —> 0, however at T y 0 it has support on the mass shell of the light field (f>. This 
feature is important, as it will be discussed in detail below, it is responsible for the dissipative processes in the 
effective action and leads to thermalization of the light field with the bath of heavy fields. It may be argued 
that an effective field theory description would only be valid for T <C M \, M 2 so that the heavy fields are not 
excited in a thermal plasma, therefore these contributions are thermally suppressed in this regime. While this is 
correct, we will see below that even when the heavy degrees of freedom are thermally suppressed, the long time 
limit will reveal full thermalization of the light field and if m</, «T«C M\ , M 2 the light field will reach a large 
population as a consequence of the decay of the heavy degrees of freedom, even when these degrees of freedom 
are thermally suppressed. The thermal suppression implies longer time scales for relaxation towards thermal 
equilibrium and the population of the particles grows during this time scale from the decay of the heavy field. 




The contribution p 2x (k 0 ,k;T) corresponds to the usual two particle cut for Q 2 > (Mi + M 2 ) 2 and is finite in 
the T —> 0 limit where it is given by (see appendix B) 


p{ko, k;T = 0) 


sign(fcp) 
87t Q 2 


Q 2 -(M 1 -M 2 ) 2 Q 2 -(M 1 +M 2 y } 0 Q 2 -{M 1 +M 2 ) 


n l 

2 r 


e[ 


; Q 2 = kl-k 2 (2.39) 


and n(ko) —> —0(—fco). The 0 function in (12.391) corresponds to the two particle threshold. This is the only 
contribution that survives at T = 0 and leads to a local and unitary effective action in this limit. The threshold 
contribution to the noise kernel has been discussed previously only at T = 0 and for Mi = Mo(42j|. 


III. T=0: LONG TIME LIMIT, LOCAL EFFECTIVE FIELD THEORY, THRESHOLDS AND 

RENORMALIZATION. 

Considering the factorized initial density matrix (12.21) and evolving it in time with the full interacting Hamiltonian 
(12.612.71) is tantamount to a sudden “switch-on” of the interaction resulting in transient non-equilibrium dynamics 
which is described by the effective action (12.311) . While the initial transients arising from the sudden “switching-on” 
of the interaction are of timely interest in a wide range of experimentally available phenomena in condensed matter 
physics, in this section we focus on the long-time behavior after the initial transient phenomena has relaxed. Our 
goal here is to study under what circumstances the long time limit after the transient dynamics has relaxed the time 
evolution of the reduced density matrix is described by an effective local Hamiltonian dynamics yielding unitary time 
evolution. We first consider T = 0, in this case the “bath” or environment is in its ground state and the spectral 
density of correlation functions of the bath degrees of freedom features only the usual multiparticle thresholds. From 
this study we draw conclusions on dissipative aspects which will highlight the influence of the medium at T y 0. 

In order to relate the long time behavior of the influence functional to the spectral representation of the bath 
correlation functions we define 

J±(fc,i) =J d 3 xe^ s ^(xy) (3.1) 

keeping the same notation for the spatial Fourier transform, secondly, for the J± with time argument t 2 we introduce 
(for the spatial Fourier transforms) 

J± (k,t)= J 


( 3 . 2 ) 
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finally we carry out the integral over £2 

r l1 . w, , , 1 ~ cos[(p 0 - fco)ii] 

/ e <Po-ko)(ti-t 2 ) dt =i 1 


(po ^ fco 


sin[(p 0 - fc 0 )fi] 
(Po - * 0 ) 


t\ —> 00 iP 


1 


Lpo - fco 


+ 7rS(po — fco) (3.3) 


where V stands for the principal part and the equivalence on the right hand side is understood when integrated with 
a smooth density of states 1 . The influence function is given by 


T[J + ,J-} = F h [J + ,J-] + T d [J+,J-], 


introducing 


A^(p 0 , k) = J ■ 


f dk ° v 

0 ,fc) 

j 2t r 

1 

O 

-se 

1 

o 

_1 


we find 


Th[J + ,J~} = -g 


d 3 k f dpo 


dt 1 e ipotl < J + (k,t 1 )J + (-k,p 0 )A > (p 0 , k) + 


and 


T d [J + , J~] 


(2tt) 3 J (2?r) J Q 

J-(k,t 1 )J-(-k,po)A < (p 0 ,k) - J + (k,ti)J~(-k,p 0 )A < (p 0 ,k) - 
J-(k,t 1 )J + (-k,p 0 )A > (p 0 ,k)\j , 

= igl ' K J /|r/ ^ i e _,po ‘ 1 |j + (fc,t i )J + (-fc,P°)p > (po,A:) 

J~(k,t 1 )J~(~k 1 p 0 )p < (p 0l k) - J + (k,t 1 )J~(-k,p 0 )p < {p 0 ,k) - 
J~(k,t 1 )J + (-k,p 0 )p > (p 0 ,k) i . 


(3.4) 


(3.5) 


(3.6) 


(3.7) 


These contributions describe very different processes: J 7 h[J + , J~] describes the contribution to the effective action 
from virtual intermediate states of the “bath” fields, as shown below, it features a non-vanishing local limit whenpoj k 
are well below the threshold for excitation of the “environmental” degrees of freedom. This contribution describes the 
“fluctuating” part of the effective action. Instead J-d[J + , J ~] describes dissipative processes. For this contribution to 
be non-vanishing, po, k must be within the region of support of the spectral density, namely above the multiparticle 
thresholds. This contribution to the effective action describes production or decay of the environmental degrees of 
freedom namely the transfer of energy and momentum between the system degrees of freedom to those of the bath, 
resulting in a dissipative contribution to the effective action. 

This interpretation becomes clear by considering the two cases (12.3712.391) for the spectral density at T = 0. 

Case a) gJ[(j>} X : hr this case for T = 0 it follows from (12.3612.371) that 


p > (k 0 , k) = — (5(fco - Wk) ; p < {k 0 ,k) = — S(ko + wu) ■ 
w k w k 


therefore 


A > (po, fc) = -—V 
2 w k 


Lpo - w k i 


A < (p 0 ,fc) = 


1 


2 w k lp 0 +w k i 


(3.8) 


(3.9) 


If the currents J ± (fc,po) only have support for k,po <C M (M is the mass of the single x held) then (13.91) can be 
replaced by —1/(2M 2 ) and 1/(2M 2 ) respectively and taken out of the integrals in (13.61) . Using (13.11) and (13.21) we find 


1 Alternatively take t\ —> oo with a convergence factor e —> 0"*" to yield 


P0~k 0 +ie ' 
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that the last two terms (13.611 cancel out, leading to 

Fh\J + ,J~] = ^ J d *xJ o dt 1 |(J + (f,t 1 )) 2 -(j-(f,t 1 )) 2 |, (3.10) 

furthermore, for po <C M it follows that p^(po, k) = 0 since the spectral density does not have support for po ^ ±Wk, 
therefore 


T d [J+,J-]= 0. 

In this case the total effective action (12.311) becomes 


(3.11) 


Seff[ft + ,ft ] = I dh I (Pxl \C Q [ft + ] + 


C 0 [ft~ 




2 M 2 


(3.12) 


This effective action describes a local effective field theory leading to unitary time evolution of the density matrix of 
the field ft with the effective local Hamiltonian 


Hef f — Hod) 


00 


9 


2 M 2 


d 3 x( J[ft}) 2 ■ 


(3.13) 


For example, with J[(j)\ = ft 2 this result simply describes an effective local “Fermi theory” by integrating out the 
heavy degree of freedom in the intermediate state, depicted in fig. ©■ The result is expected from carrying out 
the Gaussian path integrals over % without the kinetic and spatial gradient terms as befits nearly zero frequency and 
momentum transfer. It is reassuring that it emerges in the local limit of the non-equilibrium effective action. 



The term Td[J + , J~] given by (13.71) leads to non-unitary time evolution and is purely dissipative, it describes the 
coupling of the field ft to the “continuum” of the heavy field(s) as gleaned from the fact that the spectral density is 
evaluated at the “external” frequency po at momentum k. In this simple case the “continuum” is just one field mode 
and the spectral density vanishes for po,k <C M, there is no dissipation in this case and the reduced density matrix 
evolves in time with unitary time evolution, namely 

p r (t) = e""**"* p r ( 0) e iHefft . (3.14) 


Case b): for this case when T = 0 it follows from (12.3612.391) that 

p > (k o ,k;T = 0) = Q(k o )\p(ko,k-,T = 0)\ ; p < (k 0 ,k;T = 0) = 0(—fc 0 ) \p{ko,k]T = 0)| . (3.15) 


This spectral density describes a two particle continuum above the thresholds |fco| > qth = y/k 2 + + M 2 ) 2 . 

Therefore we find 


/•OO 

A > (p 0 ;k) = dko V 

J qth 


\p(k 0 ,k;T= 0) 

POO 

: A < (p 0 ; k) = / dk 0 V 

\p{k 0 ,k;T= 0) 

Po ~ k 0 

J qth 

Po + k 0 


(3.16) 


We note that as k 0 —> 00 it follows from (12.391) that \p(ko,k;T = 0)| —> l/8n therefore A ^(p 0 ,k) become log¬ 
arithmically divergent in the ultraviolet. This is precisely the zero temperature renormalization of the effective 
“current-current” vertex from the two particle intermediate state. 

As we are considering a “low energy” effective quantum field theory of a light field after integrating out the heavy 
fields, all the correlation functions of the light field are restricted to feature transferred momenta well below the two 
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particle continuum, in particular |po|, q <C qth therefore p>(po, k) = 0, and we can take po, k —)► 0 in (13.161) . to extract 
the local limit. Because A > (0;0) = — A < (0;0) we find again that the last two terms in (13.61) cancel out and using 
(liUlliHl) we find 

F„[J + ,J~] = g 2 Z g J d 3 xj^ dh |(J+(f,t 1 )) 2 - (J~(xM )) 2 
T d [J+,J-] = 0, (3.17) 

where 

Z g = f ^ |p(fc o ,0;T = 0)| (3.18) 

Jqth 

is ultraviolet logarithmically divergent. For J[<j>] = cj ) 2 the local effective vertex is shown in fig.®, and the ultraviolet 
logarithmic divergence corresponds to coupling renormalization. 

For the case J[(j)\ = <j> the local effective vertex is shown in fig.®, the loop corresponds to a self energy correction 
and includes mass and wavefunction renormalizations. To lowest order in po, k, the effective action becomes 


Seff[<i> + , 



C 0 [<t> + ]+g 2 Z g (J[<p+)) 


2 


co[<r]+g 2 z g (j[r}) 



(3.19) 


It describes unitary time evolution as (13.141) in the asymptotically long time limit after the transient dynamics has 
relaxed. 




\ / 
v * 


X2 

FIG. 3: Effective vertex in case b) for J[r/>] = ij> 2 . 


Although we focused on the limits po, k ~ 0 to extract the local limit, for po <C qth we can go beyond and expand in 
Po, k thereby generating a derivative expansion of the effective action in powers of (<9 M 4 > /(M 1 + M 2 )) that describes 
the asymptotic long-wavelength and low frequency dynamics with a local effective action. 

In the case of J{4>\ = <fi this will also contain the (finite) wave function renormalization. This is the dynamics 
expected from a local effective field theory, it is manifestly unitary because the dissipative term vanishes identically. 

This local effective and unitary description is valid for long wavelength and low frequency phenomena so that po, k 
are well below the threshold to excite the environmental degrees of freedom. The description in terms of a local 
effective field theory emerges in the asymptotic long time limit well after transient phenomena has relaxed. 

When po becomes larger than multiparticle threshold the dissipative contribution to the influence action, Tb ^ 0 
as a consequence of particle production and the time evolution of the reduced density matrix is no longer unitary. 

For T jtz 0 new contributions to the two particle spectral density arise as shown explicitly in eqn. (IB34IB35IB36I) in 
appendix ([Bj. These have support below the T = 0 two particle threshold and yield new phenomena. In particular, 
if (Mi — M 2) 2 > (pq — k 2 ) the contribution from pr> (which vanishes for T = 0) in the spectral density (12.381) (see 
IB34HB371 for the explicit contributions) will lead to dissipative processes and non-unitary time evolution which is 
studied below. 
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IV. STOCHASTIC DESCRIPTION: BROWNIAN FIELDS AND LANGEVIN EQUATION. 

In this section we study new dynamical phenomena associated with the T ^ 0 contributions to the spectral density 
with new multiparticle cuts below the two particle continuum. 

In order to establish a clear relation to a stochastic description we study the case 

J[&(x)\ = ; o(\) = xi X2, (4.i) 

which is replaced in the general effective functional eqn. (12.241) ( ( 0(x )) =0). 

Correlation functions of the fields (j>± along the forward and backward branches can be obtained by introducing 
sources hZ(x,t) linearly coupled to these fields and taking functional derivatives with respects to these sources. 
Therefore the effective action (12.31[) is generalized to 

dt j d 3 xi^C Q [4> + ] + h + (x,t)cj) + (x,t) — C 0 [4>~] — h~(x,t)(/)~(x,t)^j + T[J[(t> + ], J[4>~]]. (4.2) 

We have taken the upper time limit t —> oo because we can obtain the correlation functions of the fields at arbitrary 
times by taking functional derivatives with respect to the sources at different times. 

Therefore we consider the generating functional 

z\h±] = J Dtf p r ^f,h, oo, / 1 ± ), (4.3) 

where p r is given by eqn. (12.151) with the addition of the sources in the effective action. 

We introduce the center of mass and relative variables 

= ^((j) + (x,t)+(f>~(x,t)) ; R(x,t) = ((/) + {x,t) - (p~{x,t)) (4.4) 

and similarly for the sources 

H=^(h + +h~) ; h={h + -h~) (4.5) 

and the Wigner transform of the initial density matrix /^(O) of the (f> field 



W(4> i; IT) = j DR l e- i f d3xn ^ Ri Wpw(** + y) • (4.6) 

The Wigner transform leads to a quasiclassical description of the dynamics, it is the closest “proxy” to a (semi) 
classical phase space distribution^;, 25J, its time evolution is the Fokker-Planck equation[25|. Accordingly, the center 
of mass combination is the closest to a semiclassical description of the field[22||, in fact 

Trip p r = TrT p r ; TV R p r = 0 . (4.7) 

The inverse transform of the Wigner function is given by 

PwiVi + Y’**- T )= / £,n ie i/d3xIIi( " )ii<(S) >V(^;n i ). (4.8) 

The boundary conditions on the ^> ± path integrals given by (12.1212.131) translate into the following boundary conditions 
on the center of mass and relative variables 

^(x, t = 0) = ’Iq ; R(x,t = 0) = Ri (4.9) 

furthermore, the trace in Z entails that cj>f = (p'f leading to the constraint 


R(x , t = oo) = 0 . 


(4-10) 
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In terms of the spatial Fourier transforms of the center of mass and relative variables introduced above (14.411 . 
integrating by parts and accounting for the boundary conditions (14.91) , the non-equilibrium effective action (14.21) with 
J[4>] = (j) in eqn. (I2.24[) becomes: 

iS eff [*,K\ = dt^2{-iR_ % (%(t) + fig - ^g(t)) +^k h -k} 


> /* oo 

dt / dt f 
Jo 


: Rdt')AT k (t - f) RAt) + RAt) iEjftt - t') *At') 


d 3 xRi(x)4>(x,t = 0) 


(4.11) 


where fl 2 = k 2 Tm^,, and the last term arises after integration by parts in time, using the boundary conditions (14.91) 
and (14.101) . The kernels in the above effective Lagrangian are given by (see eqns. (12.25112.251) 1 

A4 (i-i') = 9 —[G > {kA-t') + G < {kA-t')\ (4.12) 

iT,%(t-t') = g 2 [G > (kA-t')-G < (k;t-t')\Q(t-t') = i^ k < y t-t , )Q{t~t') (4.13) 

where G < ’ > (k; t — t') are the spatial Fourier transforms of the correlation functions in (I2.25M2.251) . 

The term quadratic in the relative variable R can be written in terms of a stochastic noise as 

ex p{-^ J dt J dt' R_j:(t)Nk{t - t')Rj:{t') } = j P£exp{-i J dt J dt' £g(£).A4 _1 (£ - t')£_ u (t r ) 

+ i J dt €_z(t)Rj;{t)} (4.14) 

We now set the sources hA = 0 ( H = h = 0) to simplify the discussion, they can be included a posteriori to obtain 
the correlation functions. 

The non-equilibrium generating functional can now be written in the following form 

Z = J D^i J Dili J V't/VRV^ w{^ l ;n i )DR l e i d d 3 xR d^( n ^)-' ! ’^ t = 0 )) p[£\ x 


exp ^ — * J dt R_j:(t) 


%£) + ft^(£)+/ dt 1 E?(t-£')**;(*')-&(*) 
Jo 


where 


P[£]=expj J dt dt' CfcW-A/fc *(*-*')£_£(*')} ■ 


(4.15) 


(4.16) 


The functional integral over Ri can now be done, resulting in a functional delta function, that fixes the boundary 
condition 4^(x, t = 0) = II,; (x). 

Finally the path integral over the relative variable can be performed, leading to a functional delta function and the 
final form of the generating functional given by 


Z= U x W(4h; Tli)V*V£ V[£\6 


%(*) + + / dt' £ fc (t - £')%(*') - &(t) 


(4.17) 


with the boundary conditions on the path integral on T given by 

4/(x, t = 0) = ^j(x) ; V(x,t = 0) = ni(f) (4-18) 

where we have used the definition of T, R (t — t') in terms of £&(£ — t') given in equation (14.131) . 

The meaning of the above generating functional is the following: in order to obtain correlation functions of the 
center of mass Wigner variable *F we must first find the solution of the classical stochastic Langevin equation of motion 


+ / dt' E fc (£-f , )'Ffc(0 = %(£) 

Jo 


*£(t = 0) = *<,* ^(£ = 0 ) = n 


i,k 


(4.19) 
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for arbitrary noise term £ and then average the products of ’h over the stochastic noise with the Gaussian probability 
distribution V[^\ given by (14.161) . and finally average over the initial configurations 4' i (*);II i (£) weighted by the 
Wigner function W^i, Ilj), which plays the role of an initial (semiclassical) phase space distribution function. 

There are two different averages: 

• The average over the stochastic noise term, which up to this order is Gaussian. We denote the average of a 
functional ,A[£] over the noise with the probability distribution function P[£] given by eqn. (14.161) as 




fvzpmm 

fVCPld] 


(4.20) 


Since the noise probability distribution function is Gaussian the only necessary correlation functions for the 
noise are given by 


«&(*)» = 0 , «&(*)£_*(*')» = A4 it - t') S 3 (k - k') (4.21) 

and the higher order correlation functions are obtained from Wick’s theorem. In general the noise kernel 
Afk{t — t') ^ S(t — t') namely colored noise. From the results in appendix (1X1) we find 

iS fc (t-0 = g 2 J ^ p(k 0l k)e~ ik (4.22) 

Af k (t-f) = y J ^p(fc 0 ,/c)coth[-^]e" lfeo(t “ t ' ) (4.23) 

therefore the self energy iE and noise Af kernels obey the generalized fluctuation dissipation relation derived in 
appendix ©• 

In the literature it is usually assumed that the noise term A4(f — t') has very short range time correlations, 
namely 


ATk{t — t') °c 5(t — t') (4.24) 

this assumption entails that 

p(ko, k) coth[—~ constant, (4-25) 

in the classical limit coth[ko/2T] ~ 2T/ko this implies an ohmic spectral density, namely 

p{k 0 ,k)ock 0 (4.26) 

however, this spectral density is not compatible with that arising in the relativistic field theory from the inter¬ 
mediate state with heavy particles not even as an approximation within a range of frequencies. In particular 
the presence of thresholds implies typically power law long time tails and at finite temperature the “anomalous” 
thresholds and the non-vanishing support of the spectral densities below the T = 0 multiparticle thresholds 
clearly suggest that the ohmic spectral density cannot reliably describe the temporal correlations of a realistic 
bath. 

While such an approximation may lead to an agreement with the long time dynamics in special cases, it is 
generally a crude approximation that merits scrutiny in each case. 

• The average over the initial conditions with the Wigner distribution function yV’(’Fj, IT) which we denote as 


-4['F i ,rb] 


fDVifDIU W(^ i; n i ) 


(4.27) 
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For example for a Gaussian initial density matrix with vacuum free field correlations and anticipating a mass 
renormalization, the Wigner distribution function yields the following averages: 


rdf E = 


+ 'h,t 'h; 


-i,k i-k 2QK *> fc i ~ k 


_ n? 

n r n, r = -j- 


i.k i.—k 


n,. j: n . 


n r + ^.rii. r = n t r + ^. E n. 

IK l, — K l.K l, — K l,K l, — K l,K l, — K ' 


where are the renormalized frequencies. Renormalization is discussed below. 
The average in the time evolved full density matrix is therefore defined by 


(4.28) 

(4.29) 

(4.30) 


((-4)) == 


/2?GP[£] r 4[jqlU] 

fvam 


Calling the solution of (14.191) ^(t; £; IR), the two point correlation function of T is given by 


((*-&)*&)))■ 

It is straightforward to confirm that the equal time correlation function of the field <f> is given by 


{httW-ktt)) = ((^-k(t)^k(t))) ■ 


(4.31) 


(4.32) 


(4.33) 


where the average on the left hand side of this equation is in the reduced density matrix. 

This result is remarkably similar to the Martin-Siggia-Rose|43| formulation for stochastic classical systems. 

The solution of the Langevin equation (14.191) is obtained by a Laplace transform as befits an initial value problem. 
Defining the Laplace transforms 


pOO 

%( S ) ee / dte^^t) 
J 0 

n OO 

Cfe(s) = / dte~ st ^(t) 
Jo 


(4.34) 

(4.35) 


S (k, s ) 


dte~ st T, k (t) 


r°° p(kp, k) 

2ir J _ 00 k 0 - is 


(4.36) 


in terms of which we find 


n i,fc + s ^,fc+£fc(' s ) 

2 + n 2 + s (*, s ) 


where we have used the initial conditions (14.181) . Introducing the function G k (t) that obeys the following equation of 
motion and initial conditions 


G k (t) + n 2 k G k (t)+ [ dt' E fc (t-t , )Gfc(t , ) = 0 ; G k (t = 0) = 0; G k (t = 0) = 1 (4.38) 

Jo 

whose Laplace transform is given by 


G k (s) 


1 

s 2 + n 2 k + s (k, s ) 


(4.39) 


we find that the solution of the Langevin equation (14.191) in real time is given by 
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tEh; — He g Gk{t) + ITjg Gfc(t) + / Gk(t — t') £g(t')dt ', 




(4.40) 


the last term highlights the stochastic nature of the effective field theory. 
The real time solution for Gk ( t ) is found by the inverse Laplace transform 


G k (t)= [ 

Jc 


ds 


C 2m s 2 + ^2 + E ( fc) s ) 


(4.41) 


where C stands for the Bromwich contour, parallel to the imaginary axis in the complex s plane to the right of all 
the singularities of Gk(s ) and along the semicircle at infinity for Res < 0. The singularities of Gfc(s) in the physical 
sheet are isolated single particle poles and multiparticle cuts along the imaginary axis. Decaying states correspond to 
complex poles with a negative real part. Thus the contour runs parallel to the imaginary axis with a small positive 
real part, namely s = ito + e; — oo < w < oo with e —> 0 + . Therefore 


G k (t) = -j 


du 


2n 


(w — ie ) 2 — $1% — £#(fc, oj) — *£j(/c, w) 
which is recognized as the Fourier transform of the retarded propagator, where 


(4.42) 


£ fl (fc,w) = 


r_ 

27T 

„2 


dknV 


p{k 0 ,k) 


' - kn J 


£ i(k,u) = ■— p{u,k). 


(4.43) 

(4.44) 


where V[- • • ] is the principal part and Zi R {k, w), £/(fc, w) are even and odd functions of a; respectively as a consequence 
of the property (1A8I) which was used to write (14.4314.441) . 

In the non-interacting theory the integrand in (14.421) features poles at = ±fb, including self-energy corrections 
the positions of the complex poles are determined by 


~ - s fl(*5 w p) - i^i(k, U} p ) = 0 . 

To leading order in g 2 the real part of the poles are given by the renormalized frequencies 


w* = = klk 


Eji(k, flk) 
2 f Ik 


(4.45) 


(4.46) 


Near these poles and for small coupling, the propagator may be approximated by writing lo = uj p + (tv — w p ) and 
keeping the linear term in (w — oj p ) leading to the Breit-Wigner form of the propagator in (14.421) 




2 wb 


— i- 


(4.47) 


where the wave function renormalizaton and decay width are given by 

„2 r 


Z~ x = 1 — _ i _j_ 


9 


T k = Z 


2f Ik 47rf Ik 

Sj(fc,n fc ) g 2 Z 


dknV 


p(ko,k) 


l{n k -k 0 ) 2 } 


fb 


2 fli 


■ p(tl k ,k). 


(4.48) 

(4.49) 


In (|4~48|) = dY, R /dbj. Because £#,£/ are even and odd in oj respectively, it follows that Z,T are the same for 

both poles. The Fourier transform in (14.421) is dominated by the complex poles of the Breit-Wigner propagator, we 
find 


G k (t) 



sin[ftfc t] • 


(4.50) 
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The Breit-Wigner approximation neglects the perturbatively small contribution from the continuum along the 
branch cuts, describing the propagator in terms of simple resonances with complex poles near the real axis. The 
contribution from the continuum gives rise to power law long time tails dominated by the behavior of the spectral 
density near the thresholds, with Landau damping being the dominant contribution asymptotically for long timefl3|. 

Introducing this solution into (14.401) it is straightforward to obtain the averages of ’Fg(f) and 4/£(i)>F_j;(t) with the 
initial Wigner distribution function and the stochastic noise as in eqn. (14.311) . We find 


n 


«**(*)» = S cossin [Sl«t] 


(4.51) 


At T = 0 the spectral density of the “bath” is given by (12.391) and Q 2 = R -C (Mi + M 2 ) 2 where is the 
renormalized mass of the light <j> field, in this case T*, = 0 and only wave function and mass renormalization remain. 
After performing the average over the noise the remaining convolution integrals become similar to (13.31) in the long 
time limit. Since flj} -C (Mi + M 2 ) and only the two particle cut with threshold at Mi + M 2 remains in the spectral 
density at T = 0, in this case (T = 0 ; T^ = 0) we find that one power of Z from the terms with 4/^, H are cancelled 
by the stochastic contribution and the long time limit yields 


«*s(t)*_ s (t)» 


z 


2^L 


1 + 0 


( 9 n k \ 

\m 1 + m 2 ) 


2-i 


(4.52) 


In other words, the long time limit corresponds to the equal time two point correlation function of the field multiplied 
by wave function renormalization and with a renormalized frequency. Thus we recover the case of the local effective 
field theory limit at T = 0 emerging in the long-time limit after the transient dynamics has subsided, in complete 
agreement with the discussion in section dm 

However, for T*, 0 which is the case for T 0, dissipative effects change dramatically the long time limit. 

For Tfct 1 the contributions from the initial conditions vanish and only the stochastic part, namely the last term 
in (14.401) yields a non-vanishing contribution, it is given by 


(lilY 

V2 ng) 


dko 

(2tt) 


p(ko,k) coth[|jT 


(H«-fc 0 ) 2 +(^)' 


(4.53) 


where the dots stand for perturbatively small non-resonant contributions. 

In the narrow width approximation Tfe/H^ -C 1 and using (14.491) a straightforward calculation yields 

Z 


«%(*)*-£(*))) = 2^5 [l + 2n(fifc)] 


(4.54) 


This is a remarkable result: up to the wave function renormalization factor Z this is the equal time two point 
correlation function of a field in thermal equilibrium at temperature T and for T > it agrees with the classical 
result [44j. This result is general regardless of the origin of the decay width T. 

At T -=2 0 there are two new contributions to the spectral density (12.381) : the Landau damping term with a branch 
cut below the light cone Q 2 < 0, and the term pn with a branch cut for 0 < Q 2 < (Mi — M 2 ) 2 , their explicit 
expressions are given in appendix ([B|). The Landau damping term does not have support on the (renormalized) mass 
shell of the light field, as discussed in ref. 13] this branch cut yields a long time power law tail and is perturbatively 
small in this case. The term pr> yields a new phenomenon entirely. For M\ > M 2 m ( j >R the spectral density has 
support on the mass shell of the light particle, therefore it determines the width Tfc. This contribution to the spectral 
density describes the decay process xi ~(► X 2 <t> of the heaviest particle in the medium and the inverse process. These 
processes lead to a damping of the fluctuations of the <f> and the thermalization of the light particle <fi with the heavy 
particles in the bath. Thus the coupling of a light degree of freedom to a bath of heavy degrees of freedom that 
can decay into the light field turns the <j> particles into quasiparticles [38|, [45| with a finite lifetime that describes the 
thermalization of the latter. The decay of the heaviest environmental field builds the population of the light particle 
over a time scale t t h — 1/T until it thermalizes with the bath. 


A. Quantum kinetic interpretation of the <j> width: 

The physical reason by which the decay xi X'Z'P aR d the inverse process leads to thermalization of the light <j> 
particle with the bath of heavy degrees of freedom can be understood from the following quantum kinetic analysis. 
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The interaction Hamiltonian 


Hi = 9 J d 3 xxi{x)x2{x)cj)(x) . 


with Mi > M 2 ?7i0 describes the decay and inverse processes in the medium 

Xi X 2 <t>- 


(4.55) 


(4.56) 


The transition amplitude for the gain term of <fi particles from the decay of the heavier particle in the medium, namely 
Xi X2<P is given by 


M fi 


( 2 ti -) 6(w^ - w™-n q ) 


gain 


W 


2 uP 2wl 2) 2FL 


«i(l + ^2)(1 + N^q)) 


k ~~ k+q 

where n\^ are the equilibrium distribution functions for the heavy fields xi ,2 

m = n(w J^) ; n 2 = u(wi 2) ) , p=\q+k\ 


(4.57) 


(4.58) 


and Nif, is the out of equilibrium distribution for the light particles </>. The total gain probability per unit time is given 

by 


d-Bgain _ fj~ 


dt 


(27r) 


d 3 k X 




2Q q x J ( 2 tt ) 3 


2toP 2 w _ 

K k-\-q 


J 2 ) - n i(l + 7*2X1 + N^q)) 


(4.59) 


The loss term from the inverse process X 2 4> —■> Xi is obtained from the above expression by the replacement 

Tii —> (1 + 7i 1 ) ; 1 + n 2 —> n 2 ; (1 + N^^q)) —> N^^q) (4.60) 

leading to the quantum kinetic equation for the population of </> particles 

dX! ( p(qj dPg a [ n (//’loss 


namely 


dt 


dN^q) _ g 2 tt f d 3 k n *) f 


dt 


dt ’ 


dt 4f l q J (2n) 3 

These processes are depicted in fig. ©■ 


( 1 ) ( 2 ) 
W k W k+q 


7*1 (1 + 77 2 )(1 + N^q)) - (1 + ni)n 2 N c/> (q) 


(4.61) 


(4.62) 
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FIG. 4: Quantum kinetic interpretation of the decay width: the decay process yi —>• X 2 </> and its inverse lead to a build up of 
the population of <j> particles in the medium. 


This kinetic equation may be written as 

dN^q.t) 

dt 


= (1 + N4q,t))rf-N^q,t)rl 


(4.63) 


Inspection of (IB12I) along with the relation (IB 141) shows that r> is given by the second term in (IB12I) with qo = Q(g) 
and T< is given by the third term of p < (q0 = fl q ,q). Furthermore it is straightforward to confirm that 


r> — p PClqT'< 
q 


(4.64) 






























leading to the final form 
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where 


dN c j > (q' 1 1) 
cit 


= -7(9) N^t) - N eq (q ) 


N eg (q) 


1 

e pa. { q) _ y 


and 


7(9) = r> - r 


< 

Q 



p D (qo =n q ,q) 


(4.65) 


(4.66) 


(4.67) 


where pd(<Zo = q) is the second term of the spectral density of the bath (12.381) evaluated qo = f l q . The solution of 

USD is 


N<t,{q;t)) = N eq (q) + 


N^q; 0 ) 


N eq {q) 


£- 7 ( 9 )* 


(4.68) 


Since pn this is the only term in the spectral density with support at qo = it follows that 

2 

7 ( 9 ) = Jy p(qo = n q ,q) = r q (l + 0{g 2 f) (4.69) 

where in the last equality we used the results (14.4414.491) . Consequently, 

T<=T q n(n q ) ; T> = T q [l + n(fi,)] . (4.70) 

This is one of the important results of this article: the decay of heavy fields in the medium into the light field 
leads to the thermalization of the light field with the heavy degrees of freedom. The dynamics of thermalization is 
non-unitary and is manifest in the dissipative kernels which are determined by the support of the spectral density of 
the bath on the mass shell of the light degree of freedom. This results in that the light field is described as a resonance 
whose width is precisely the rate of approach to thermal equilibrium with the bath of heavy particles. 

The equivalence between the results from the stochastic description and the quantum kinetic equation confirms the 
arguments of ref. [46] within a different framework and approach. 


B. J[(f>] ^ (f>: multiplicative noise 


The case J[4>] = 4> studied above highlights the stochastic nature of the effective action in a clear manner and 
exhibits the generalized fluctuation-dissipation relation between the self-energy of the light field and the correlation 
function of the stochastic noise. However the effective Langevin description is particular to the linear coupling to the 
bath. 

For a generic coupling we now implement a Kramers-Moyal expansion in the relative coordinate I 22 I [24l. l25j . writing 
in terms of the Wigner center of mass (\t f ) and relative coordinates ( R ) as in (14.41) . 

J[^] = J[4/] ± | J'[T] + ^ J"[4/] + • • • (4.71) 

Inserting this expansion into the influence action (12.331) and into the effective action (12.311) we find 

/»00 

iS e ff[^,R) = dt 
Jo 

dt J d 3 x J dt' J d 3 a/ji R(x, t) J'[^t(x, t)] M{x — x'\ t — t') R(x', t') J'[4'(x / , t')] 

+ R(x,t)J'['$/(x,t)]i'E R (x~x';t — t')J['I>(x',t')]'j + f d 3 xRi{x)^/(x, t = 0) 



d 3 ‘. 


"j — iR(x,t) — V 2 4/(x, t) — j 


(4.72) 
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Again, the term quadratic in R can be written in terms of a stochastic noise as 

expj — i J dt J d 3 x J dt'J d 3 x'R(x,t)J'['I>(x,t)]Af(x — x';t — t')R(x',t')J'['I>(x',t')]^ 
J 2 ?£expj — i J dt J d 3 x J dt' j d 3 x' ^{x,t)N~ l {x — x\t — t')^(x',t') + 
i J dt J d 3 x £(x, t)R(x, t) J'\$?{x, f)]j. 


We find that the exponential term in (14.1611 is now given by 


exp|—i J dt J d 3 x R(x,t) t) — t) + m^(x, t)+ 

J dt! j d 3 x' J'\$>{x, t)]H R (x — x']t — t')J\$>{x' ,t')\ — J'\$>(x, t)]£(:r, t) 


where 


(4.73) 


(4.74) 


V [5] = exp | — — J dt J d 3 x J dt' J d 3 x’ £(x, t)Af~ 1 (x — x', t — t')£(x l , f')^ . (4.75) 

The Langevin equation of motion for T that follows from the effective action is now 

^{x, t) — V 2 5'(x, t) + m^(x, t) + J'[ty{x, t)] J dt' J d 3 x T, R (x — x 1 -, t — t') J\$!(x!, t 7 )] = J'\^>{x, t)] £(x, t) , (4.76) 

with 


«£(£>*)» = 0 ; ({£(x,t) £(x,t'))) =N(x-x l ,t-t'). (4.77) 

Therefore we see that when J[(j>] ^ </> we obtain also a stochastic description, but with a multiplicative Gaussian noise. 

From the time evolution of the Wigner function, or generating functional (14.16|) we could implement the steps 
detailed in refs.[20i. [22l - l25l | to obtain a Fokker-Planck equation, however the generating functional (14.161) in terms of 
the influence action when generalized to include external sources as in eqn. (14.31) contains already all the information 
necessary to obtain correlation functions at equal or different times by functional derivatives with respect to the 
sources. This is akin to the Martin-Siggia-Rose framework of stochastic classical phenomena [43j| . 


V. THE QUANTUM MASTER EQUATION: 

In the Master equation approach[l4|, [l5( the time evolution of the density matrix in considered in the interaction 
picture. With the full density matrix p(t) given by eqn. (12.51) it follows that 

p I {t)=e iH ° t p{i)e- iH o* ; H 0 = H 0 [<j>\ + H 0 [ X ] (5.1) 


whose time evolution obeys 


pi{t) = -i[H I (t),pi(t)\ (5.2) 

where Hj(t) is the interaction Hamiltonian given by eqn. (12.71) in the interaction picture, namely the fields 
(j>(x,t);x(x,t) evolve in time as free fields. The solution of (15.11) is 


pi(t) = p{0 ) - i 


dt' pi(t')] 


This solution is inserted back into (15.21) leading to the iterative equation 


(5.3) 


pi(t) =-i[H I (t),p(0)\ - / [i7/(f), [#/(*'): P/(*')]] dt 1 , 

Jo 


(5.4) 
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the next step relies on a series of assumptions , the first being : Factorization : 


Pi(t) = pitit) ® p x (0), (5.5) 

namely that the heavy degrees of freedom (the “bath”) remain in thermal equilibrium. Taking the trace over the 
X degrees of freedom an normal-ordering the interaction following eqn. (12.321) the first term on the right hand side 
of eqn. (15.41) vanishes upon taking the trace over the environmental degrees of freedom, and we find the evolution 
equation for the reduced density matrix for <f> in the interaction picture, 


Pi<t>(t) 



d 3 x' 


J(x) J(x') pic/,(t') G > (x 


x') + pi<j,{t') J{x') J(x) G < (x — x') 


- J(x) Pi4>{t') J(x')G < (x - x') - J(x')pi^{t') J(x) G > (x - 



(5.6) 


where again J( x) = J[(f>(x)\, and we use the shorthand convention x = (x,t) ; x' = (, x',t' ). The hat on J is to 
emphasize that these are (composite) operators in the interaction picture, G > ,G < are given by (12.2512.261) and we 
suppressed the subscript c (connected) because the operator O has been normal ordered. At this stage a Markov 
approximation is usually invoked by replacing pi^it') —► pi^ft) taking it outside the time integral, this approximation 
is justified in weak coupling. For example consider the first term in (15.61) . it can be written as 


2 J(x,t ) f —j^-pi(j,(t')dt' ; 'H(t')= f J(x',t") G > (x — x',t — t")dt h 
Jo dt J q 


which upon integration by parts yields 


(5.7) 


- g 2 J{x,t)'H{f)p I(t> {t) + g 2 J(x,t) J 'H(t') dt' (5.8) 

in the second term dpi^t') / dt' oc g 2 so this term yields a contribution that is formally of order g A and can be neglected 
to second order. The same analysis can be applied to all the other terms in (15.61) with the conclusion that in weak 
coupling and to leading order ( g 2 ) the Markovian approximation pi<j>(t') —» pi<j>(t ) is justified. 

Therefore in the Markov approximation the quantum master equation becomes 


pi4,(t) 



d 3 x 


d 3 x' 


J(x) J{x') pi<j>{t) G > (x — x') + pi<t,{t) J{x') J(x) G K (x — x') 


- J{x) pi<t,(t) J(x')G < (x - x') - J(x') pi<j>{t) J(x)G > (x 



(5.9) 


The relation between this master equation in operator form and the effective action of the previous section which 
is cast in terms of functional integrals in the field basis is not a priori evident. However we can infer the equivalence 
from the result given by eqn. (12.341) . The first step is to recall that the variables with + correspond to the forward 
time evolution branch, namely from the time-ordered evolution, whereas those with — correspond to the backward 
time evolution branch, namely, anti-time ordered evolution. Because the time evolution of the density matrix is 
determined by /7(t)p(0)t/ _1 (t) with U(t) the forward time evolution operator (time ordered) and C/ -1 (t) the backward 
time evolution operator (anti-time ordered), operator insertions with + go before the density matrix, and operator 
insertions with — go after the density matrix. Therefore when multiplying the density matrix the equivalence becomes 

A+B+ -> ABp ; A~B~ -> pAB ; A+B~ -> ApB ; A~B + -> BpA (5.10) 

Therefore this equivalence applied to (15.91) states that the first term in (15.91) is identified with J + (x)J + (x')G > (x — 
x'), the second with J~(x)J~(x')G < (x — x'), the third with — J + (x)J~{x')G < {x — x') and the fourth with 
— J~(x)J + (x')G > (x — x'). From these identifications we immediately recognize that these are precisely the terms 
in the time derivative of the influence function (12.341) . which describes the interaction of the light particle with the 
bath of heavy fields. This “dictionary” establishes the direct correspondence between the influence function and the 
master equation approaches: the solution of the quantum master equation 15.611 in the field basis and in the Markov 
approximation is the influence action \2.33\) . 

The factorization assumption (15.51) is equivalent to obtaining the influence action from the correlation functions 
of the bath in the initial density matrix (of the bath) and exponentiating the result in a cumulant expansion, as 
manifestly exhibited in the result (12.21112.291) for the influence action. 
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A. Local effective field theory: effective Hamiltonian. 

In order to understand how a local effective field theory emerges from the quantum master equation, and to establish 
contact with the results of section let us consider the interaction gJ[4>] X with only one heavy field (\) of mass 
M mp at T = 0. Writing the correlation functions G & as in eqn. (IA1IA2I) with the spectral densities given by 
(13.81) . defining 


J d 3 xJ[x,t]e ikS = j dp 0 e- ipot J{k,p 0 ), (5.11) 

and carrying out the integral in t' in (15.91) in the long time limit t —> oo with a convergence factor e —> 0 + , we find 


Pi<t>(t) = -ig 2 


d 3 k 


(27 t ) 3 2 w k 


dpo / dq 0 e 


-i(qo+Po)t j J{k } p 0 )J(-k,q 0 ) p I(j> (t) p I(j) (t) J(k,q 0 )J(-k,p 0 ) 


(qo ~Wk + ie) 


(qo + w k + ie) 


J(k,po) J{-k,q 0 ) J(k,q 0 ) pip(t) J(-k,p 0 ) 


(qo + w k + ie) 


(qo ~w k + ie) 


(5.12) 


If the frequency qo and momentum k transfer are qo,k -C M we can set qo ~ 0 in the denominators and replace 
Wk = M. In this low energy limit the last two terms in (15.121) cancel each other out, and from (15.111) we obtain 


Pip(t) = ~i [Hl ff (t),p I4 ,(t)} , 


(5.13) 


where 




2 M 2 


d 3 xJ 2 [(f>(x,t)] . 


(5.14) 


Therefore we recover the local “current-current” limit with Hamiltonian evolution for the reduced density matrix in 
the interaction picture in complete agreement with the result (13.131) from the influence action approach. This result 
confirms that obtained from the influence action in section dm In particular, the cancellation of the last two terms 
in (15.121) in the “local limit” is the same as the cancellation of the last two terms in (13.61) . furthermore in this local 
limit the limit e —> 0 can be taken safely without yielding a delta function, such term would yield a purely real, namely 
dissipative, contribution to the quantum master equation. The cancellation of the last two terms and the vanishing of 
the dissipative contribution in agreement with the results of section m are a direct manifestation of the equivalence 
between the influence action and the quantum master equation discussed above. 


B. J[<P] = </> Lindblad Master equation. 


The equivalence with the influence function and the results of the stochastic Langevin equation in particular the 
dynamics of damping and thermalization become more clear by studying the case J[<fi\ = <j> to compare with the 
results of section m■ In the interaction picture and in terms of the spatial Fourier transform in a volume V 


<t>& t) = -4y 

Q 


(5.15) 


where 


Mt) = 


y/XT q L 




(5.16) 


and the operators a^\ a[_ - do not depend on time. At this point we invoke yet another approximation: the “rotating 
wave approximation”: in writing the products <j>(x,t) c/)(x',t') in (15.61) there are two types of terms with very 
different time evolution. Terms of the form 


at a 9 - 




(5.17) 
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and its hermitian conjugate are “slow”, and terms of the form 


4°'-f 


,2 in 


ft t ) 


axa,- 


_^ — 2i€l q t e — iSlq(t — t') 


(5.18) 


are fast, the extra rapidly varying phases lead to rapid dephasing and do not yield resonant (energy conserving) 
contributions. These terms only give perturbatively small transient contributions and are discussed below. Keeping 
only the slow terms which dominate the long time dynamics and neglecting the fast oscillatory terms defines the 
“rotating wave approximation” ubiquitous in quantum optics. We will adopt this approximation and comment later 
on the corrections associated with keeping the fast terms. Implementing the Markov approximation pi^lt') —> pj^it), 
and the rotating wave approximation (keeping only terms of the form ad a, a cd) using the spectral representation 
of the correlators (12.3511 with the property p < {— ko, k) = p > (k 0 , k) and carrying out the spatial and temporal integrals 
we obtain the Lindblad form of the quantum master equation, 


where 






,Pi<fi{t) 


r£(t) 


r£(t) 


4 h pi <*(*) + Mt) a l a k - 2a k Mt) 4 


k fe 


;4 + Mt) H a \ - 2 4 Mt) a k 


Rk(t) = 


47rf lb 


dk 0 p(k 0 ,k ) 


1 — cos[(flfc — k 0 )t] 


(f^fc ko) 


(5.19) 


(5.20) 


211 


dk 0 p{k 0 , k ) [1 + n(k 0 )] —^41 

k J 7r(i2fe - k 0 ) 


(5.21) 


r<u\ 9 2 f ai (i n n ^m[(Q k - k 0 )t] 

r * {t) = 2n k ] dfcop ( fco ’ fc)n ( fco) n(n k -ko) 


(5.22) 


The second and third lines in (15.191) are called the dissipator , these are non-Hamiltonian, purely dissipative terms 
that cannot be written in terms of a local hermitian effective field theory. 

Taking the long time limit at this stage using the results on the right hand side of eqn. m , would lead to 


R(t) t —> oo 9 [ dk 0 V 

47T ilk J 


P(k 0 ,k) 

k Q 


Y, R (n kl k) 

2 f Ik. 


= SCl k , 


(5.23) 


which is recognized as the renormalization of the frequency from eqns. (14.4314.461) . Similarly in the long time limit 

„2 


r fc(0 t^-oo p(flfe, k) [l + n(flfe)] = r 


(5.24) 


Tfc(i) OO ^-p(n k ,k)n{n k ) = r< 


(5.25) 


r k( t )~ r k( t ) t^oo ^-p(n k ,k) = T k (5.26) 

where r^,Tfe are precisely the rates in the quantum kinetic equation (14.6314.701) and coincide to leading order in 
g 2 with the result in eqn. (14.491) . Therefore in the long time limit the first term in the quantum master equation 
(15.191) is identified as a local contribution to the effective Hamiltonian (in quantum optics this term is referred to as 
the Lamb shift), and the second and third contributions correspond to the dissipative terms in the Lindblad master 
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equation. In particular at T = 0 the long time limit leading to (15.2415.251) would vanish because at T = 0 the 
spectral density does not h ave support on the mass shell of the light field, and features a two particle threshold at 
|fco| = yjk 2 + (Mi + M 2) 2 » flfc- Therefore we will proceed to analyze the consequences of the Lindblad master 
equation (15.191) without taking the long time limit at this stage. 

For any interaction picture operator A associated with the light field </> 

j t (A) = Tv{A pi*® + A , (5.27) 

where the average ((• • •)} = Tr(- • • )pi^,{t). Since ag , at are time independent in the interaction picture, we find 


where 


^ / t \ 

S ( “E> 


— iR(t) — 


r k (ty 


(°k> 


iR{t) — 


r k {t) 


K-> 


r fe (t) = r>(i)-r<(i). 


Using the asymptotic integrals 


t 



dko 


p(kp,k) 
k kp) 


sin(Ufe — fco) t 

(f h - k 0 ) t 


00 



- p{k 0 ,k ) ~ 
.(U fc - ko)- 


(5.28) 


(5.29) 


(5.30) 


f°° ,, P(kp,k) 

7-oo ° ( k 0 -n k ) 2 

we find in the long time limit 


1 — cos [(fc 0 ~ Ufc)t] 


t —> 00 irtp(£l k ,k) + / dui' V 


p(k 0 ,k) 


H^k-k 0 y 


rt—too 

/ R(t')dt' = t 
Jo 


(5.31) 


(5.32) 


r k (t') . r k 1 

dt' = — t + z~ x -1 


(5.33) 


where Z“ 1 is given by eqn. (14.48[) . The emergence of the wave function renormalization is a consequence of keeping 
the time dependence in T?(t) in (15.191) until the full solution is obtained. Taking the long-time limit too hastily in 
the Lindblad equation does not allow to extract the wave function renormalization asymptotically. 

Therefore in the long time limit we find 


(«fe)(i) 


Z(a s )( 0 ) e~ i5Qk 



<a£)(t)~Z<4>( 


(5.34) 


Introducing these solutions into {4> k {t)) by replacing (a^); (at) into (15.151) and recognizing that + 5Q k = we 
find in the long time limit 


(<%(*)> 


Ze~ 


(^(°))cos[f2ft] + 


Ml 


sin [£l k t\ 


K-(°)> = !<«*)> 


(5.35) 


this long time solution is the same as that obtained with the influence functional given by eqn. (14.51[) . 

To establish the relation with the influence action result (14.541) we now consider the evolution equation for averages 
of bilinears: 


*Ml = -T k {t)N k (t) + r<(t) ; N k (t) = (a\a % ), 


(5.36) 
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d 

dt 


a_g) = - 2 iR k (t) -T k {t) (a £ a_g) 


^<4 a 4) = [ 2iR k(t) - r fc (t)] (at a f _-), 


(5.37) 


Taking the long time limit in Tj7(f) —> T^ - the rate equation (15.361) agrees with eqn. (14.651) from quantum kinetics, 
neglecting the (perturbative) contribution from wave function renormalization, the solution of (15.361) in the long time 
limit is given by 

N k (t) = N eq , k + [JV fc (0) - N eq , k ] e~ Tkt (5.38) 

which describes thermalization when M\ > M 2 m^>, namely when (at least) one of the heavy fields in the medium 
can decay into the light field, in agreement with the influence action, stochastic and quantum kinetics approaches to 
leading order in g 2 . Furthermore, the solutions to the equations (15.371) for the bilinears become in the long time limit 


{a k a_ k )(t) = (c 


a _fc)(0) e 2iRki e rict , 


(5.39) 


and its hermitian conjugate, again neglecting wave function renormalization perturbative contributions. Combining 
these results we find in the long time limit 


(^feW^-fcW) 1 + 2 N eq ,k 


(5.40) 


which again exhibits a stationary thermal correlation function in agreement with the results of the influence action to 
leading order in g 2 . The renormalization of the frequency Cl k —> fi R can be included by writing the free field theory 
of the light field (j> in terms of the renormalized mass and introduce a counterterm in the perturbation, in the long 
time limit the Hamiltonian term R k (t) in the Lindblad equation will cancel the frequency counterterm. 

The solutions ( 15 . 381 ) and ( 15 . 401 ) emerges in the long time limit provided the spectral density has support on the 
mass shell of the light field (f> namely if the decay process \i ~> X 2 ^ is available. At T = 0 for example the spectral 
density has a two particle threshold above the 4> mass shell and T k = 0. In this case there are no resonant terms, 
if the long time limit is taken (hastily) in Tjf’(f) (| 5 . 2115 . 221 ) then the dissipative contribution to the Lindblad master 
equation would vanish identically. 

However, keeping the time dependence of T^T ( t ) reveals important aspects that are missed if the long time limit is 
taken at the outset, these emerge from the solution of the rate equation (15.361) with the time dependent rates. It is 
given by 


N(t) = 


IV(0) 


r fc (*') efo 


r k (t")dt" dt 1 


(5.41) 


The long time asymptotic behavior is obtained from the identities (15.331) . if T k ^ 0 the secular term dominates at 
long times and memory of the initial conditions is lost after a time 1 /T k leading to the result in eqn. (15.381) . 

On the other hand, in the case when T k = 0 as in T = 0, the long time limit of the time integrals are non-secular 
and given by (14.481) . where Z~ x — 1 oc g 2 , therefore to leading order in g 2 we can neglect the exponential terms in 
(15.411) and assuming that the field 4> is initially in the ground state, namely IV(0) =0 we find that asymptotically 

N(t 00 ) oc Z~ l — 1 oc g 2 . (5-42) 


This result has a clear interpretation: after the transient dynamics subsides the coupling to the bath of heavy particles 
has produced virtual excitations of the 4> held, and the “true” (or dressed) ground state is a linear superposition of the 
“bare” ground state and multiparticle excitations. The wave function renormalization measures the overlap between 
the bare and the dressed states and the number of “bare” particles in the “dressed” state. 

An important conclusion from the equivalence with the stochastic description is the identification between the 
“rotating wave approximation” and the Breit-Wigner approximation for the propagator (14.471) which replaces the full 
propagator by a complex pole which describes reliably the long time limit. 

Counterrotating terms: 

In the derivation of the Lindblad master equation (15.191) we neglected terms of the form 

a ^ a _ e -2in fe t . t „t 2iQ k t -iQ k (t-t') 

k -k 0 ’ k -k 


(5.43) 
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The time integral over t' can be carried out as in the Lindblad form yielding contributions of the form 
a k a -k k)pi^{t) etc. The contribution of these terms to the equations of motion for linear or bilin¬ 

ear forms of a, a' are straightforward to obtain, they do not yield terms that grow secularly in time because the rapid 
dephasing of the oscillatory terms average out in the time integrals. A simple analysis shows that these terms yield 
perturbatively small subleading contributions of the form Sflk/^k ; T k/^k as compared to those obtained from the 
Lindblad form which captures the secular growth in time because of the resonances and describes the leading behavior 
in the long time dynamics. 

In conclusion, the “counterrotating” terms always yield perturbatively small contributions that are bound in time. 
They are negligible in the case when the spectral density of the bath has support on the mass shell of the light 
field (resonant terms), and are perturbatively small in the case when it does not, in agreement with the perturbative 
corrections obtained from the influence action approach in this case, (see eqn. (14.5211 ). 


VI. DISCUSSION 

We have established the relation between the quantum master equation and the influence action and compared the 
evolution of expectation values and correlation functions from the two approaches. Each approach has advantages 
and disadvantages that merit a discussion. 

• Influence action: the effective action (14.21) and the generating functional (14.31) yield the time evolution of the 
density matrix for the system in the presence of external sources and functional derivatives with respect to these 
sources yield any arbitrary correlation function. The stochastic nature of the time evolution is more explicitly 
manifest in this approach and the correlations of the bath enter in a more explicit manner in the statistical 
and stochastic averages. An important advantage of the influence action is that it allows to obtain correlation 
functions at different times from functional derivatives with respect to the external sources. 

• Quantum master equation and Lindblad form: This approach is more useful to obtain directly the time 
evolution of slow averages, namely averages of operators that do not evolve in time in the interaction picture. 
The connection with the quantum kinetic description is more readily recognized in this approach as is manifest 
in the results (15.361 lATTl) . The general time evolution of averages of operators is given by eqn. (15.271) . the first 
term of which requires the time evolution of the density matrix, not just its time derivative. Therefore in order 
to extract the time evolution of averages of arbitrary operators it is necessary to obtain the full time evolution of 
the density matrix, only the time evolution of averages associated with operators that do not depend on time in 
the interaction picture are simple(r) to obtain. Correlation functions at different times are not straightforward 
to obtain from the quantum master equation either in Lindblad form or the more general form including 
the counterrotating terms. The relationship between the general form of the quantum master equation m 
discussed in section © and the time derivative of the influence function (12.341) clarifies that the influence action 
is the solution of the quantum master equation for the reduced density matrix in the interaction picture in the 
field basis. Furthermore this equivalence also sheds light on the nature of the factorization assumption in the 
master equation formulation. Therefore the effective action (12.311) and in particular the generating functional 
in terms of the Wigner transform (14.171) give the time evolution of the reduced density matrix. 

• Thermalization: In our study thermalization emerges as a consequence of the decay of a heavy particle in 
the bath into another heavy particle of smaller mass and a light particle. This case is actually quite general: 
in principle one would expect that the application of effective actions to low energy physics (or long time, 
long wavelength phenomena) is a result of tracing over a large number of heavy degrees of freedom of various 
masses and couplings to the light degrees of freedom. Therefore the particular case in this study is, quite 
likely, representative of the general situations in which an effective low energy description is sought. When 
the temperature of the environment is comparable to the energy scale of the heavy degrees of freedom, these 
will be excited and present with large population in the plasma. However, thermalization of the light degree 
of freedom will occur even when the temperature of the environment is much smaller than the masses of the 
heavy fields. Obviously the population of these fields will be exponentially suppressed, and this is reflected 
in the decay rate T^ through the spectral representation (see appendix [Bj, however, this entails that the time 
scale to full thermalization t ~ 1/T becomes very long, nevertheless thermalization will ensue and the light 
“system” will eventually reach a stationary, thermal state on this time scale. The non-perturbative nature of 
this result is noteworthy: if the temperature is much smaller than the masses of the heavy fields, the time scale 
to thermalization is very long and the secular growth of the fluctuations become very large. In the long time 
limit even for small coupling and large masses, the fluctuating part of two point correlation function of the light 
field at equal time approaches the thermal distribution (15.401) . As discussed in detail in the above sections, 
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this process of thermalization is necessarily described by the non-local dissipative contributions to the effective 
action. 


VII. CONCLUSIONS AND FURTHER QUESTIONS. 

In this article we studied the emergence of an effective field theory out of equilibrium from the open quantum 
system perspective where a light field <j> of mass - the “system”- is in interaction with a bath or environment 
at temperature T, taken to be either one scalar heavy field \ with mass M or two scalar heavy fields XI 1 X 2 

with a hierarchy of masses Mi > M 2 m^. We obtain the reduced density matrix of the field <f> by tracing out the 
environmental degrees of freedom up to second order for interactions of the form </[</>] O[xj] where J, O are in general 
polynomials of the fields. The time evolution of the reduced density matrix is determined by the influence action 
which is generally non-local and describes dissipative processes and non-unitary time evolution of the reduced density 
matrix. For the case when the interaction is of the form g J[<fi\ \ we show how the familiar “current-current” effective 
field theory emerges from the influence function in the low frequency, long wavelength limit at long times after the 
transient dynamics subsides. 

When the environment contains two heavy fields XuX2 the spectrum of environmental fluctuations is very rich 
and leads to a wealth of dynamical phenomena. For an interaction g<px 1 X 2 a Wigner (semiclassical) transform of 
the effective action yields a stochastic description of the dynamics described by a Langevin equation with non-local 
dissipative and additive Gaussian noise kernels that obey a generalized fluctuation dissipation relation. When the 
interaction is non-linear in <f> the noise is multiplicative. At T = 0 the spectral density of the bath only features a 
two particle threshold well above the mass shell of the light field and the influence action of the bath yield mass and 
wave function renormalizations, the asymptotic long time behavior corresponds to a renormalized local effective field 
theory. 

For T / 0 new “anomalous” thresholds emerge and the spectral density of the bath features support on the mass 
shell of the light field that describes the decay xi —> Xify- This “in medium” process is manifest as non-local dissipative 
contributions to the effective action which ultimately lead to the thermalization of the light field with the bath. Even 
for T < Mi, M 2 when the population of heavy degrees of freedom is thermally suppressed, the light field thermalizes 
at long times. We show that the dissipative contributions and the thermalization of the light field from the decay 
of the heavy fields in the bath are in complete agreement with quantum kinetic results. We argue that this result is 
quite general in that an environment with heavy degrees of freedom will include a hierarchy of very massive states 
which upon interacting with the light fields will result in their decay into the light degrees of freedom. This decay 
necessarily implies that the spectral density of the bath will have support on the mass shell of the light particle, which 
in turn leads to a dissipative contribution to the effective action and eventual thermalization of the light field. 

We also obtain the quantum master equation up to second order in the interaction and show directly that its solution 
in the field basis is precisely the influence action and elucidate the nature of the various approximations invoked in 
the quantum master equation approach. For the case g<fixiX 2 we obtain the Lindblad form of the quantum master 
equation under precise approximations which are also understood from the equivalence with the influence action. 
It contains a Hamiltonian and a dissipative term whose coefficients depend explicitly on time. We show that this 
time dependence is crucial to obtain consistently important renormalization aspects. The asymptotic time evolution 
obtained from the Lindblad master equation with the time dependent coefficients is in complete agreement with that 
obtained from the influence action, the stochastic Langevin equation and the results from quantum kinetics. 

While the influence function, Langevin and quantum master equation approaches all agree in the long time limit, 
they offer advantages and disadvantages which we discussed in detail. In particular, the influence action approach 
leads directly to the stochastic Langevin description, while the stochastic nature is not a priori evident in the quantum 
master equation. Furthermore, the infuence action approach when augmented to include external sources allows to 
obtain correlation functions at different times which is more difficult in the master equation approach. 

These results offer a note of caution on the application of effective field theories in a finite temperature (and likely 
a finite density) environment such as the early Universe, as dissipative effects arising from the influence of heavy 
environmental fields lead to non-local and non-unitary stochastic dynamics of the light degrees of freedom. 

Further questions: 

The time evolution of the initially prepared density matrix leads to the entanglement of the system and environment 
degrees of freedom. After tracing over the environment degrees of freedom, the influence action describes the time 
evolution of the reduced density matrix and is a direct result of the entanglement between the system and the 
environment. The Von-Neumann entropy associated with the reduced density matrix is the entanglement entropy 
(although alternative definitions of the entanglement entropy are available 17|). Several questions remain to be 
explored: a) how to extract the entanglement entropy from the influence action, b) whether the time evolution 
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of the entanglement entropy has a quantum kinetic interpretation, as for example an //-theorem, c) what are the 
manifestations of the stochastic nature of the effective action upon the entanglement entropy. 

Our ultimate goal is to study the emergence of effective field theories in cosmology, in particular during the 
inflationary stage, under the assumption that there are heavy degrees of freedom with mass or energy scales larger 
than the Hubble scale that are traced over leading to an effective description in terms of a single “light” scalar field. 
In inflationary cosmology there are novel processes associated with the lack of a global time-like Killing vector [HU] 
that are not available in Minkowsky space time and preclude a spectral representation of correlation functions of the 
environmental degrees of freedom. In particular there are no kinematic thresholds and quanta of a field can decay in 
quanta of the same held (5^. f53j|. or of a heavier held suggesting that the effective held theory description obtained 
upon tracing over degrees of freedom with mass scales much larger than the Hubble scale during inhation will always 
be dissipative and non-unitary. The results of this study will be reported elsewhere. 
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Appendix A: General fluctuation dissipation relation: 


Because the initial density matrix p x (0) = e P H ° X is translational invariant and 0(x, t) = e lH ° xt 0(x, 0) e we 
can write 


G > (x — x']t — t') = (0(x, t)0{x', t')) 
G < (x — x';t — t') = {0(x', t')0(x, t)) 


d 4 k 

d A k 
(27r) 4 


p>(k, fc 0 )e _ifeo(t “ t ' ) e^- 5 ") 
p<(k, fc 0 )e" ifco(t “ t ' ) . 


Introducing a complete set of eigenstates of H 0x it follows that 


p > {k 0 ,k) 

p < (k 0 ,k) 


J2 e pEn {n\O^)\m){m\O_^)\n)S(k 0 - (E n - E m )) 
J2 e ~ l3En (n\0_j:(0)\m){m\Oi:(0)\n)5(ko- {E m - E n )) 

rX \ ) rn . ri 


(Al) 

(A2) 


(A3) 

(A4) 


where 0%( 0) ; 0_%( 0) is the spatial Fourier transform of the (composite) operators 0(x,t = 0) ; 0(x',t' = 0) 
respec tively. Upon relabelling m f->- n in the sum in the definition (IA4I) we hnd the Kubo-Martin-Schwinger relation fl2L 


p < (k 0 ,k) = p > (-k 0 ,k) = e /3k °p > (k 0 ,k ) (A5) 

where we have used parity and rotational invariance in the second line above to assume that the spectral functions 
only depend of the absolute value of the momentum. The spectral density is defined as 

p(k 0 ,k) = p > (ko, k) - p < (k 0 ,k) = p > {k 0 , k) [l - e~ pk °] (A6) 

therefore 

p > (k 0 ,k) = p{k 0 ,k) [l +n(k 0 )} ; p < {k 0 ,k) = p{k 0 ,k) n(k 0 ). (A7) 

Furthermore, from the first equality in (IA5I) it follows that 

p(-k 0 ,k) =-p(k 0 ,k). (A8) 


In terms of the spectral densities we find 

G > (x — x 1 -, t — t') — G < (x — x r ;t — t') 


d 4 k 

Wf 


p{k 0 , k)e 


—iko(t—t) ^ik-(x-x') 


(A9) 
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where 


G > (x — x;t — t') + G < (x — x';t — t') 



/c(m - 

JC(k 0 ,fc)e- ifeo(t - t ' ) e i£ -^- £ ' ) 


(A10) 


/C(fc 0 ,fc) 


p(fc 0 ,fc)coth[-^] . 


(All) 


This is the general form of the fluctuation dissipation relation, with 1C proportional to the kernel for the noise 
correlation function (see eon. (14. 121) 1. Note that p(ko,k ) is odd whereas lC(ko,k) is even in fco- 


Appendix B: Spectral densities: 

1. Case a): single x field of mass M 

The finite temperature correlation function for a single \ field is 

G>{2-i-,t-0 = <*(*, tMi'X)) = /j Ii±^>l + ntC 

where Wk = y/k 2 + M 2 . Relabelling k — > —k in the second term (1B1I) can be written as 

d A k 


e iw k (t-t') e ~ik-(x-x') 


G > {x — x'-,t — t') = 




p > (k 0 ,k) e 


—iko(t—t') gik-(x-x') 


where 


p > (k 0 ,k) =— S(k 0 - w k ) - 5(k 0 + w k ) [l + 7i(fc 0 )l . 
w k L 1 


(Bl) 

(B2) 

(B3) 


In a similar manner we find 

G < (x-x'-,t-t') = {xix 1 ,t')x{x,t)) = 

which is written as 


f d k j [l ~i~ Tl(wk)~\ iwk(t—t ') -ik-(x—x') _j_ n {.Wk) ik-(x—x r ) 


(27r) 3 I 2wife 


2 Wk 


G < (x — x'\t— t') = 


d A k 

W? 


p < (k 0l k) e 


—iko(t—t') gik-(x-x') 


where 


The spectral density 


p > (k 0 ,k) = — S(ko - Wk) - S(ko + Wk) n(ko). 


w k l 


p(k 0 , k) = p > {ko, k) - p < {ko , k) = — 6(k 0 - Wk) - 6(k 0 + Wk) 


w k l 


(B4) 

(B5) 

(B6) 

(B7) 


2. Case b): Xi; X 2 > Ml > M 2 
Defining in four vector notation K = (fco, k) 

(x j y) Xj (y)) = J^i^(K)e- iK ^- y) ; 3 = 1,2 


(B8) 
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with 


A j(k 0 ,k) = 


(1 + n(w'j?’)) 5(k 0 - w^) + n(w { k j) ) S(k 0 + w ( k JI ) 


,.(?)> 


i 


U) = Jk 2 + M 2 


(B9) 


With O(x) = Xi( x )X 2 (x) it follows that 


G>{x-y) = ( 0 (x) 0 (y)} = ( X x(x)xi(y))(X2(x)x2(y)) = j / ^A l( tf) A a (P) 


^0p > (q 


from which we obtain 


p > (qo,q) = 


d 3 kdk 0 


Ai (k 0 ,k) A 2 {qo - k 0 , \q- k |). 


(BIO) 


(Bll) 


Using m we find 

P>{q ° ,q) = 2/(2?^ 


d 3 k 1 


(1 + ni)(l + n 2 ) 5(q 0 - - w^) 


+ (1 + ni)n 2 S(q 0 - w < j}' 1 + w^) + (1 + n 2 )n 1 6(q 0 + - w + nin 2 S(q 0 + + uj£ 2) ) j (B12) 


ni=n(tu[. } ) ; n 2 = n(w^ 2) ) ; p=\k + q\. 


where (after relabelling k —> —k) 

..1 - ; n 2 = «Ui( 2 ) 

From the general relation 

p < {qo:q) = p > (-qo,q ) 

obtained above (see eqn. (1A5|I ) and p(qo, q) = P > (qo , q) — P < (qoi q )> we finally find 

1 


(B13) 

(B14) 


7 t f d 3 k 
p{qo,q) = 77 


(1 + m + n 2 ) 5(q 0 - w[ 1] - w^ 2) ) - d(q 0 + w k + w 2 ) 


2 J (2tt) 3 w W w ( 2 ) 

+ (n 2 - ni) [<5(g 0 - + w4, 2) ) - S(q 0 + w ( k 1} - w {2) ) 


(B15) 


In the first line in (12.381) the terms with n 2 can be simplified by relabelling k —> — k — q in the k integral, with 

—> Wp ^ ; w p 2 '* —> w^ k \ then the term proportional to n 2 in the first line becomes of the same form as that for 

ni but with the replacement M\ o M 2 . In the second line a similar relabelling makes the term proportional to n 2 
similar to that of n\ upon M\ 0 M 2 . Therefore 

p(qo,q) = p {0 \qo,q) + (p (/) 0zo,«) +M X M 2 ) + (^p {II) (q 0 ,q) + M 1 o M 2 ) (B16) 

where 

P (0) (<7o,<7) = ^2 si g n (<?o) f ( f) fc (2) <HM ~ - w {2) ), (B17) 

J W^ Wp 

P (I) (qo,q ) = Y^ si g n (®) f (f) fc (2) S(\qo\- w { k ] - wf>), (B18) 

P (H) (qo,q) = Y^sign(go) f -^-^y n^) <5(|go| + w { k ] - w {2 >), (B19) 
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Obviously, p^°\qo,q represents the zero temperature contribution. Let ui = and z = with 

Q-\~K 


Then, we have 


d 3 k = k 2 dkd(cos(0)) ; 


d(cos(9)) dz 

w (2) - k q 

q+k 


P ( °\qo,q) 

p {I) (qo,q) 

p (II) (qo,q) 


l 


167T 2 

1 

167T 2 

1 

167T 2 


/*oo /*~ 

sign (go) / dcj / S(\q 0 \- uj - z)dz 

J Mi Jz~ 

nOO r Z + 

sign(go) / n(cj) dco / 5(\qo\ — co — z) dz 

J Mi Jz~ 

s-OO nZ + 

sign(g 0 ) / n(u) dto / 5(\q 0 \ + u - z) dz , 

J Mi Jz~ 


(B20) 


(B21) 

(B22) 

(B23) 


where 


z ± = ifij^qf^M 2 = \Juj 2 + q 2 ± 2 q^/to 2 - Mj 2 - (M 2 - M 2 ). 


(B24) 


Without loss of generality we assume that Mi > M 2 . 

For the integrals in (IB21IB22I) to be non-vanishing, it must be that 


z < |< 70 1 - < z + . 


(B25) 


A simple analysis shows that for large q 0 the curve |go| — intersects both z ± (k), as |g 0 | diminishes, both 
intersections occur with z~(k) until they coalesce, leading to the condition Q 2 = qfi — q 2 > ( Mi + M 2 ) 2 . We find 


P (0) (go, q) = S1 ^ 2 0) B(Q 2 - (Mi + M 2 ) 2 ) j [Q 2 - (Mi - M 2 ) 2 ] [Q 2 - (Mi + M 2 ) 2 


(B26) 


where 


p {I) (qo,q) = 


sign (gp) 
8nfdq 


0(Q 2 - (Mi + M 2 ) 2 ) In 


1 - e~P w + 
1 - e~d w - 


w±(qo,q) = ^{|go|« ± g^W} 


(B27) 

(B28) 

(B29) 


a = Q 2 + Ml - M| ; a 2 - AQ 2 Ml = Q 2 - (Mi - M 2 ) 2 Q 2 - (Mi + M 2 ) 2 . 

The analysis for (IB23II follows the same steps, now the integral with the delta function is non-vanishing for 

z~(k) < |g 0 1 + < z + (k) . (B30) 

There are two different cases: i): |go| < g (or Q 2 < 0) or ii): |g 0 | > g. In case i): there is intersection only with z~(k) 
and the range of integration in to is 


where 


£ < w < 00 , 


£(go,g) = 2 |^ 2 | I \qo\a + q\Ja 2 + d\Q 2 \M 2 


(B31) 


(B32) 
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In case ii): there are two intersections for 0 < Q 2 < (Mi — M 2) 2 and no intersections for Q 2 > (Mi — M 2 ) 2 , combining 
the two cases we find 




We summarize this result as 


1 


Li — 


-0(Q 2 )0((M 1 -M 2 )"-Q 2 ) In 


1 - e~ 0w + 
1 - e~P w - 


where 




8ir(3q 1 Ll — J 


+ Mi±* M 2 


(B33) 


p(q 0 ,q;T) = p LD (qo, q; T) LD 0(-Q 2 ) + p D (q 0 , q\ T) 0((Mi - M 2 ) 2 - Q 2 ) 0(Q 2 ) 

+ p 2x (q 0 , q ; T) 0(Q 2 - (Mi + M 2 ) 2 ) ; Q 2 = q 2 0 - q 2 (B34) 


(B35) 


/ rp' sign(go) f 

= In 


1 - e~P w + 

1 - e~P w - 


+ Mi <h> M 2 


(B36) 


P 2 X {qo,q]T) = 


sign(g 0 ) / 
87rQ 2 I 

sign (go) f , 


Q 2 - (Mi - M 2 f || Q 2 - (Mi + M 2 ) 2 
1 - e -P w + 


1 - e~P w - 


+ Mi o M 2 


(B37) 
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